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ABSTRACT 

We present the PLANET photometric dataset* for the binary-lens 
microlensing event MACHO 97-BLG-28 consisting of 696 I and V"-band 
measurements, and analyze it to determine the radial surface brightness profile 
of the Galactic bulge source star. The microlensed source, demonstrated to 
be a K giant by our independent spectroscopy, crossed an isolated cusp of the 
central caustic of the lensing binary, generating a sharp peak in the light curve 
that was well-resolved by dense (3 — 30 minute) and continuous monitoring 
from PLANET sites in Chile, South Africa, and Australia. This is the first time 
that such a cusp crossing has been observed. Analysis of the PLANET dataset 
has produced a measurement of the square root limb-darkening coefficients 
of the source star in the / and V bands; the resulting stellar profiles are in 
excellent agreement with those predicted by stellar atmospheric models for K 
giants. The limb-darkening coefficients presented here are the first derived from 
microlensing. They are also among the first for normal giants by any technique 
and the first for any star as distant as the Galactic bulge. Modeling of our light 
curve for MACHO 97-BLG-28 indicates that the lensing binary has a mass ratio 
q = 0.23 and an (instantaneous) separation in units of the angular Einstein 
ring radius of d = 0.69. For a lens in the Galactic bulge, this corresponds to 
a typical stellar binary with a projected separation between 1 and 2 AU. If 
the lens lies closer (i.e., in the Galactic disk), the separation is smaller, and 
one or both of the lens objects is in the brown dwarf regime. Assuming that 
the source is a bulge K2 giant at 8 kpc, the relative lens-source proper motion 
is = 19.4 ± 2.6 kms -1 kpc -1 , consistent with a disk or bulge lens. If the 
non-lensed blended light is due to a single star, it is likely to be a young white 
dwarf in the bulge, consistent with the blended light coming from the lens itself. 

Subject headings: binaries — stars: atmospheres and fundamental parameters 
— Galaxy: stellar content — techniques: microlensing 



*Based on observations at: Canopus Observatory, Tasmania, Australia; the European 
Southern Observatory, La Silla, Chile; and the South African Astronomical 
Observatory, Sutherland, South Africa. 
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1. 



Introduction 



The PLANET (Probing Lensing Anomalies NETwork) collaboration uses a 
longitudinally-distributed network of telescopes in the southern hemisphere to perform 
densely-sampled, precise photometric monitoring of Galactic microlensing events with the 
express goal of searching for anomalies, especially those that may betray the presence 
of extrasolar planets ( [Albrow et al. 1998|) orbiting the lenses. Multiple lenses, such as 
lensing binary stars or planetary systems, generate measurably anomalous light curves if 
the background source comes close to a caustic generated by the lens; direct transits of a 
caustic can cause sudden and dramatic changes (amplification ^ 10 on time scales of hours) 
in the apparent brightness of the source. Three such microlensing events were monitored by 
PLANET in its 1997 Galactic bulge observing season; here we describe our results for one 



Although typical subgiants and giants in the Galactic bulge have radii that subtend 
only ~1— 10 microarcseconds (//as), the dramatic change in magnification in the vicinity 
of caustics can selectively magnify the spatial structure of the source star by huge factors, 
momentarily acting as a very high-resolution, large-aperture telescope trained on the 
background source. These relatively rare but spectacular lensing events are particularly 
important because they can yield otherwise unobtainable information about the lens and 
source. 

Generally speaking, the only physical parameter that can be deduced from non- 
anomalous microlensing light curves is the characteristic time scale of the event, defined 
as the time required for the lens to move an angular distance across the observer's sight line 
to the source equal to the angular Einstein radius 9e 



where M is the total mass of the lens and D^, D$, and D^s denote the observer- lens, 
observer-source and lens-source linear distances, respectively. 

In caustic crossing events, photometric monitoring can allow more information to 
be obtained, including the relative proper motion fi of the lens-source system and — if 
the coverage is complete over the caustic — the surface brightness profile of the source 
star. During a caustic transit, the radial surface brightness profile, or limb- darkening, 
of the source influences the shape of the resulting light curve ( [Witt 19951) . Photometric 
data of sufficient precision and temporal resolution can thus be inverted to deduce the 
limb-darkening profile of source stars as distant as the Galactic bulge ( [Bogdanov fc| 
Cherepashchuk 1996| ; |Sasselov 1997] ), including possible chromatic effects induced by 



of these, MACHO 97-BLG-28. 
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differential limb-darkening at different wavelengths ( |Bogdanov fc Cherepashchuk 1995a 



1995b; |Valls-Gabaud 1995| , 1998; |Gould fc Welch 1996| ). Time-resolved spectroscopic 



monitoring with very large aperture telescopes may reveal even more information, including 
detailed information about stellar atmospheric physics ( Loeb fc Sassciov 1995 ; |Gaudi fc 



Gould 1999| ; |Heyrovsky, Sasselov fc Loeb 1999"|) and metallicity ( [Lennon et al. 1996|) of the 



sources. 

One of the most important observational constraints on models of stellar atmospheres 
has come from limb-darkening measurements of the Sun, a relatively cool dwarf. Constraints 
on the atmospheres of other types of stars are more difficult to obtain. Measurements of 
radial profiles, or more specifically, limb-darkening parameters, have been attempted for 
only a handful of stars (see [Scholz 1997| and references therein). Eclipsing binaries, lunar 



occultations and interferometry have been used to resolve and measure the surface structure 
of nearby stars. 

The classical method of using the light curves of eclipsing binaries has provided 
accurate limb- darkening measurements or consistency checks in only a few cases (Popper 
1984; Andersen 1991; Ribas et al. 1987), but at least one suggests that a non-linear 
limb-darkening law is indicated (Milone, Stagg & Kurucz 1992). Direct HST imaging of the 
bright, nearby M2 supergiant Betelgeuse has resulted in estimates for its limb-darkening 
and other surface irregularities (|GilIiland fc Dupree 1996| ). Interferometric techniques also 
require very large, bright stars (Hanbury Brown et al. 1974; Dyck et al. 1996; Mozurkewich 
et al. 1991; Mourard et al. 1997), where surface structure is often due to spotting or 
dust envelopes rather than limb-darkening (Roddier & Roddier 1985; Wilson et al. 1992). 
Although the interferometric signal is often incompatible with a uniform stellar disk, with 
the exception of extreme supergiants like Betelgeuse ( |Cheng, Hege fc Hubbard 1986| ; |Burns 



et al. 19§7|), theoretical limb- darkening laws are adopted in order to arrive at better stellar 



radius determinations (Quirrenbach et al. 1996; Hajian et al. 1998). Nevertheless, the 
technique holds promise for limb-darkening determinations as it begins to probe into the 
normal K giant regime ( Pyck, van Belle fc Thompson 1998Q . Lunar occultations have 
been used to resolve stellar surface structure (Bogdanov & Cherepashchuk 1984, 1990, 
1991; Richichi & Lisi 1990; Di Giacomo et al. 1991), but the dust associated with these 
supergiants complicates limb-darkening estimates here as well (Richichi et al. 1991, 1995) 
and surface profiles from theory are generally assumed so as to derive better diameters 
(White & Feierman 1987; Richichi et al. 1998). 

The precision of microlensing in mapping stellar surface brightness profiles can thus 
make an important contribution to stellar atmospheric theory by providing limb-darkening 
and spectral information on stars too faint or too small to be studied by traditional 
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methods. We report here our precise and temporally dense photometric observations of the 
MACHO-detected microlensing event 97-BLG-28, which when combined with our spectral 
typing and modeling have produced the first limb-darkening measurement via microlensing, 
and the first resolution of surface structure in a star so distant (~8 kpc). 



2. PLANET Observations of MACHO 97-BLG-28 

The Galactic bulge microlensing event MACHO 97-BLG-28 was alerted by the MACHO 
Collaboration on 29 May 19970 as a normal event expected to reach peak brightness 
near 10 June 1997. PLANET began to observe the event directly after the initial alert as 
part of its normal monitoring program. Within hours of the anomalously rapid increase 
in its brightness on 14 June 1997, PLANET^ and MACHO/GMAN independently issued 
electronic alerts announcing 97-BLG-28 as a likely caustic-crossing binary event. Dense 
photometric PLANET monitoring continued over the peak of the light curve and for about 
6 weeks thereafter, with additional data taken near the baseline magnitude in October 1997 
and at the beginning of the 1998 season. A pre-caustic spectrum was taken on 31 May 1997. 



2.1. Photometric Monitoring 

Three PLANET telescopes participated in the monitoring: the Dutch/ESO 0.92 m 
at La Silla, Chile, the SAAO lm at Sutherland, South Africa, and the Canopus lm near 
Hobart, Tasmania. The full PLANET dataset consists of 686 data points taken in 1997 
(508 I-Band: 247 La Silla, 128 SAAO, 133 Tasmania; 178 V-Band: 98 La Silla, 38 SAAO, 
42 Tasmania) and another 10 data points (5 I-Band; 5 V-band) collected at the SAAO lm 
in March and April of 1998. For details concerning the telescopes, detectors, and field 
sizes, refer to Albrow et al. (1998). The data were reduced with DoPHOT ( |5chccht~ 



Mateo fc Saha 19931) using fixed position catalogs and four stars in the field chosen to be 



stable, moderately bright and relatively uncrowded as a relative flux standard ([Albrow et 



al. 1998|). Instrumental magnitudes were calibrated against contemporaneous observations 



of Johnson-Cousins, UBV(RI)c E-region standards (Menzies et al. 1989) made at SAAO 
at the beginning of the 1998 season. As we describe in §3.1, individual zero points for 
photometry at different PLANET sites allow us to refer all our multi-site data to the 



IOts 



MACHO alerts are posted at http://darkstar. astro. washington.edu 



11 PLANET anomaly alerts are posted at http:/ /www. astro. rug. nl/^planet 
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SAAO natural system, and thus to the standard system. Measurements taken in 1998 and 
calibrated in the manner described above produced baseline magnitudes of the source star 
of Johnson V = 17.91 and Cousins / = 15.66, accurate to 0.05 mag. 

Whenever the event was fainter than / « 14.7, measurements derived from frames 
with poor image quality (FWHM > 2.2") were clearly and systematically fainter than 
those from moderate to high quality images, and had a scatter that was substantially 
larger than the DoPHOT-reported error bars. In order to avoid systematically biasing our 
modeling, measurements with I ;> 14.7 (corresponding to V 17) and FWHM> 2.2" were 
excluded from the fitting procedure to extract lens and source parameters. In addition, 
one particularly poor image that produced measurements that deviated by 0.75 mag for 
constant stars as bright as I = 14 was also excluded. This high-quality subset of our dataQ 
consists of 586 data points (431 I-Band: 247 La Silla, 130 SAAO, 54 Tasmania; 155 V-Band: 
98 La Silla, 41 SAAO, 16 Tasmania) and is displayed in Fig. 1. 

2.2. Spectroscopy 

On 31 May 1997, prior to the anomalous peak of MACHO 97-BLG-28, red and blue 
spectra were taken of the source star by Dr. M. Sahu using the high-throughput EFOSC 
(ESO Faint Object Spectrographic Camera) on the ESO 3.6m at La Silla, Chile. The camera 
has both imaging and spectroscopic capabilities, and can hold 5 gratings simultaneously so 
that multiple spectral resolutions and bandpasses can be easily configured. 

A direct image of the field was first obtained with an exposure of 10 sec. The 512 x 
512 pixel Tektronix CCD provided a 5.2' x 5.2' field of view. No filter was used in order to 
minimize any possible offset between the slit and the direct image. A slit of 1 .5 width was 
used for all the observations, and was placed in the correct parallactic angle. Observations 
were conducted through thin cirrus clouds in seeing of ~l"0. The blue spectrum was 
obtained at 04:40 UT using the B150 grating, which covers 3900-5300 A, with a resolution 
of about 5A. The red spectrum was obtained at 05:25 UT using the O150 grating, which 
covers the wavelength range 5200 - 6800 A, with the same resolution. The exposure time in 
both cases was 2400 sec. 

The image processing software packages MIDAS and IRAF were used to reduce the 
spectral data. After bias subtraction, the frames were flat-field corrected with an average 
of 5 flat-field images of the dome illuminated with a tungsten lamp. To ensure good sky 



12 This PLANET dataset for MACHO 97-BLG-28 is publicly available at frttp://www.astro. rug.nl/~planct 
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subtraction, sky levels were determined from both sides of the spectrum. The resulting 
one-dimensional spectrum was then wavelength calibrated with a He-Ar spectral lamp 
and flux calibrated against the standard star LTT 4364. Shown in Fig. 2 is the resulting 
geocentric (lab frame), combined spectrum for the source star of 97-BLG-28. Obvious 
cosmic rays have been removed. The spectral resolution is ~ 5A, corresponding to a 
velocity resolution of about 260 km at the strong Nal absorption feature near 5890A. 
The heliocentric correction for the velocities is 9.5 km s~^, which is negligible compared to 
the resolution of the spectrum. 

2.3. Radius and Spectral Type of the Source 

Comparison of the flux-calibrated spectrum shown in Fig. 2 with spectral synthesis 
codes (FITGRID and SPECFIT tasks of IRAF) yields the best match to a K2III star of 
solar metallicity, effective temperature T = 4350 K and surface gravity log g = 1.9. The 
strong Nal line is not typical of a star of this type, but the equivalent width of the line is 
entirely consistent with interstellar absorption along this long sight line to the bulge source 
star at celestial (J2000) coordinates a = 18 h 00 m 33.8 s and 5 = -28°01'10", corresponding 
to Galactic I = 2.46° and b = —2.36°. Some ambiguity remains in the spectral solution 
depending on the extinction to the source, but the best fit for the reddening of the spectrum 
and the radius of the source are consistent with that derived from an independent analysis 
of the color- magnitude diagram of field, described below. 

The radius of the source star was determined by considering its position in the 
(V — I)o — Jo color-magnitude diagram (CMD) relative to red clump (RG) stars in the 
Galactic bulge. Our CMD was derived from observations taken at the SAAO lm, calibrated 
by reference to E-region standard star baseline data in 1998. 

The absolute / magnitude of the peak of the RC luminosity function in the bulge has 
been determined to be Irq = — 0.23± 0.03 (Stanek & Garnavich 1998; Paczyhski & Stanek 
1998). We have fitted the /-magnitude distribution of the RC stars in our sample using 
the same function as Stanek & Garnavich (1998). Based on the distribution of stars in 
CMD, we find that there is an equal probability that the 97-BLG-28 source star is in the 
RC or on the red giant branch (RGB). Possible values for the average field reddening were 
determined by fitting the observed RGB to the Bertelli et al. (1994) isochrones for ages of 
3.1 and 10 Gyr and metallicities [Fe/H] = —0.4, 0, +0.4, resulting in the range of reddening 
values 0.8 < E(V — I) < 1.2. Stanek's (1998) reddening map of the Galactic plane indicates 
a gradient of ~ 0.2 mag across the 3x3 arcmin SAAO field; interpolation at the position 
of the microlensing event gives E(B — V) = 0.88, corresponding to E(V — I) = 1.10. 



- 8- 



For each value of plausible reddening, we used the known RC absolute magnitude to 
rescale our CMD, treating the distance modulus as a free parameter. For each isochrone fit 
to the RGB for a given age and metallicity, other isochrones were also plotted to allow for 
the possibility that the source star may differ in age or metallicity from the RGB mean. In 
all cases for which a reasonable fit (A(V — I)o < 0.1 mag) was found for the source being 
on the RGB or in the RC, the source radius was calculated from the appropriate isochrone. 
Reasonable consistency emerged among the 25 resulting radius values, with all in the range 
13 < R*/R Q < 20, and a mean value of R*/R & = 15 ± 2 (1-er uncertainty). This result will 
be used together with the modeling in §|5] to estimate the relative proper motion of the lens 
with respect to the source. 

In Fig. 3, we show our CMD for the field of 97-BLG-28, dereddened with our fitted 
value of E(V — I) = 1.095 (Ai = 1.56) over the whole of the field. The position of the 
observed (blended) source star at baseline is shown as the filled circle in Fig. 3, where 
it is assumed that the source suffers the same reddening as the field mean. Our derived 
dereddened color of (V — I)q = 1.16 from CMD considerations is in excellent agreement 
with the (V — 7)o = 1.14 (Bessell 1990) expected for our K2 giant spectral typing of this 
source. 

3. An Extended-source, Binary-lens Model for MACHO 97-BLG-28 

The abrupt change of slope in the light curve near HJD - 2449719 = 896.4 evident in 
Fig. 4 suggests that the source star has crossed the caustic structure of a binary lens. The 
single nearly-symmetric finite peak near HJD - 2449719 = 895.6 further suggests that the 
extent of the source is resolved by the observations and that the source passed behind an 
isolated cusp of the caustic geometry on a nearly perpendicular trajectory. This lensing 
geometry is indeed borne out by our modeling of the light curve for 97-BLG-28, as can 
be seen in Fig. 5. To reproduce the features in the observed light curve it was necessary 
to include lens binarity, source blending, extended source size, and source limb-darkening 
in this model. Together with allowances for possible photometric offsets between the 3 
observing sites in both bands, this requires a total of 19 fit parameters, as we now explain. 

We label each subset of the data by observing site and filter. Let k denote the 
number of such data subsets. A model with a point-source and point-lens will require 
3 + k parameters: £e = #e//-£, the (Einstein) time required for the lens to move an angular 
distance 6*e (Eq. 1); uq, the smallest lens-source angular separation in units of 8e] t , the 
time at which this smallest separation occurs; and the set of k m l , the multi-site, multi-band 
baseline magnitudes. Individual baselines are required for each site because although all 
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data of a given band are referred to the same secondary field standards, we often find 
that small discrepancies remain in our multi-site photometry that appear to be related to 
differences in detector resolution and average seeing conditions (Albrow et al. 1998). 

Parametrization of a binary lens composed of two objects with mass Mi and M 2 
requires an additional 3 parameters: d, the angular binary separation in units of 9e] 
q = M 2 /Mi, the binary mass ratio; and a, the angle between the line from M 2 to M 1 and 
the direction of source proper motion (relative to the lens). These parameters and the 
direction of proper motion are chosen so that Mo>0,0<a!<27r,0<g<l, and the 
midpoint of the lens system is on the right hand side of the moving source as viewed by the 
observer. For a binary system, the smallest angular separation u between lens and source 
now refers to the midpoint of the lens system, and the angular Einstein radius #e refers to 
the total mass M = M 1 + M 2 . 

The baseline magnitudes include all light, including any that may be due to any 
unresolved, non-lensed light that may be confused with source light in the observations. 
To account for the possibility of such photometric confusion, blending parameters fj 
(j = 1 . . . n) must be added for each of the observed n wave bands to characterize the 
fraction of the total light contributed by the source alone at baseline. No blending in wave 
bands j corresponds to fj = 1. The extended size of the source is characterized by an 
additional parameter = 0*/#e, the angular radius of the source in units of the Einstein 
radius. Finally, the light intensity profile of the source may be limb-darkened rather than 
uniform over the stellar disk. We adopt a limb-darkening law of the form 

J A (#) = J A (0) [l - c x (1 - costf) -d x (l- V^)} , (2) 

where d is the angle between the normal to the stellar surface and the line of sight. The 
coefficients C\ and d\ are wavelength dependent, requiring an additional 2n parameters. 

The inclusion of all these effects thus requires 7 + k + 3n fit parameters. For the 3 
PLANET sites, each observing MACHO 97-BLG-28 in V and I, k = 6 and n = 2, thus 
necessitating 19 parameters in the full model. 



3.1. Fitting the Model to the Data 

We fit three different extended-source binary-lens models, with either a uniform source 
or a 1- or 2-parameter limb-darkened source, to the combined, high-quality dataset of 
586 CCD frames (§2.1) using the formal uncertainties reported by DoPHOT. Details of 
calculating the light curves for such models can be found in the work of Dominik (1998). The 
multi-site, multi-band baselines were allowed to vary independently in the fitting process; 
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this photometric alignment resulted in relative multi-site offsets of 0.007 — 0.065 mag in our 
best model. 

In general, it is a difficult task to find a minimum in high- dimensional parameter space. 
However, in this case, nearly optimal values for some parameters could be found by first 
searching in lower- dimensional subspaces. The six baselines could be estimated from the 
latest data points and then fitted with a point-lens model together with the parameters 
u 0i to, tE and the two blending parameters, using the data points outside the peak region. 
We thus began fitting the binary-lens extended-source models only after we had good 
guesses for 11 of the parameters. Parameters from the uniform source fit then provided us 
with good initial estimates for a total of 15 parameters, requiring only 2 or 4 additional 
parameters when including limb-darkening. 

Our best fit was achieved for a 2-parameter limb- darkening model (LD2), and yielded 
a x m in °f 1913 for the 567 degrees of freedom (d.o.f). This model is displayed in the left 
panels of Fig. 4 and its fit parameters are listed in Table |l|. If the LD2 model is indeed 
the best representation of the data, then the reduced Xmm/d-O-f- = 3.374 is an indication 
that the DoPHOT uncertainties are on average underestimated by a factor ~1.8, consistent 
with our previous experience ( [Albrow et al. 19"9"% ). 

The corresponding source trajectory and caustic structure for the LD2 model shown 
in Fig. 5 illustrate that the cusp caustic of the binary lens just sweeps over the limb of 
the extended source, the first clear observation of such a cusp crossing. Such a crossing 
differentially magnifies different portions of the source as a function of time (top panel, 
Fig. 6). Because the central ring crosses almost directly over the cusp, it experiences the 
greatest magnification. The fraction of light contributed by concentric rings (of equal area) 
on the stellar disk therefore also varies during the anomaly (bottom panel, Fig. 6), allowing 
the surface profile, and thus limb-darkening, of the source to be deduced. For the specific 
geometry of 97-BLG-28, the contrast between the fraction of light contributed by the inner 
and outermost of 10 equal-area rings is a factor ~3 during the cusp crossing. The PLANET 
dataset is well-sampled everywhere during the cusp crossing except in the 12-hour period 
near HJD - 2449719 = 895. 



4. Limb-Darkening of the Bulge K Giant Source Star 

In addition to our best fitting model with square root limb-darkening coefficients (LD2), 
Table |] also lists the parameters for two other models: a fit with a linear, 1-parameter, 
limb-darkening law (LD1), which yields xLm °f 1930 for 569 d.o.f., and a fit with a 
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uniformly bright source (UB), which yields Xmm °f 3255 for 571 d.o.f. Under the assumption 
that our 2-parameter limb-darkening model is indeed the best representation of the data, 
the large x 2 compared to the number of degrees of freedom must be attributed to a 
misestimation of the experimental uncertainties. This is not unexpected since our previous 
experience in crowded fields indicates that DoPHOT formal uncertainties are smaller than 
the true scatter in the photometry of constant stars by an amount that depends on the 
crowding of the star (Albrow et al. 1998). In order to assess conservatively the differences 
in xinn betw een models, therefore, we first rescale the formal DoPHOT uncertainties by 
^1913/567 = .374 PS 1.8, which is tantamount to the assumption that the LD2 model 
is perfect and would yield a true Xmm = 567 if the "true" photometric uncertainties were 
known. The rescaling is global, preserving the relative uncertainty between points, and 
conservative, allowing for the possibility that part of the Ax min between models can be due 
to the underestimation of photometric uncertainties by DoPHOT. The rescaled Xmm °f the 
LD1 and UB models are thus 572 and 965, respectively. 

The significance of the limb-darkened LD1 and LD2 models over the uniform brightness 
UB model is enormous: using the rescaled uncertainties, A% min = 393 between the LD1 and 
UB models, and A% min = 398 between LD2 and UB. This precipitous drop in Xmm with the 
addition of the limb-darkening parameters leaves no doubt that limb-darkening has been 
detected to an exceedingly high degree of confidence in the PLANET photometry of this 
event. The enlargements of the light curve in Fig. 4 show how precisely the 2-parameter 
limb-darkening model (LD2) reproduces the detail in the cusp and limb regions in both 
the / and V bands. On the right of the same figure, these sections of the light curve are 
overplotted with the uniform source model (UB). The multi-site photometry is aligned 
automatically by the fitting procedure; this photometric alignment is somewhat different 
for the two models. Even given the freedom to photometrically realign the data from 
different sites relative to one another, the UB model fails to fit the shape of the convexity 
at peak and the gentle slope of the curve as the limb of the star egresses. The failure of 
the uniformly bright model to reproduce the structure in the light curve during the limb 
crossing can be seen most dramatically in Fig. 7, where the residuals of the LD2 and UB 
models are shown. 

If the smaller normalized Xmm °f the LD2 model compared to the LD1 models is not 
physically significant, but instead only due to the freedom of adjusting the additional d\ 
parameters, we would expect Ax min to be distributed as x 2 with a number of degrees of 
freedom equal to the number of additional parameters (Press et al. 1986, Chapter 14). 
The significance of the LD2 model over the LD1 model is thus determined by the 8.2% 
probability of obtaining by chance a Ax min > 5 improvement with two additional degrees of 
freedom, corresponding to a marginal 1.7a detection of a surface profile that deviates from 
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the 1-parameter limb-darkening law in both bands. This marginal improvement indicates 
that the inclusion of additional profile fitting parameters beyond those included in the LD2 
model is unlikely to result in a significantly better fit; the LD2 model contains the all the 
information that we are able to pull from this dataset about the source profile. 

The large difference in between the limb-darkened and uniform bright models is 

not due solely to the cusp-crossing portions of the light curve. Since the model parameters 
are correlated, even those parameters unrelated to the source profile differ between the 
uniform source model and the limb-darkened models, as inspection of Table [I] indicates. 
Interestingly, the limb-darkened models always required a smaller photometric offset 
between sites in both bands. The uniform brightness model apparently attempted to match 
the observed shape of the light curve at the limb by adjusting the photometric offsets 
slightly, hindering the y 2 performance of the UB model elsewhere in the light curve. The 
clear signature of limb-darkening is revealed only with high precision data during the few 
hours when the limb of the star is grazing the caustic; high-quality data over the whole of 
the light curve is required, however, to obtain a full and accurate microlensing solution. 



4.1. Robustness of the Measurement 

Another test of the reliability of our limb-darkening parameters is provided by fits 
we performed on the full dataset, including the poorer quality frames, but now estimating 
the photometric uncertainties empirically so that the size of the error bar scales with 
the overall frame quality. Specifically, we set the uncertainty of the event magnitude on 
a given frame equal to the average deviation of all similarly-bright constant stars from 
their average magnitude. We also eliminated the baseline points taken in 1998 in order to 
test whether the availability of a long temporal baseline was important to our conclusions 
about limb-darkening. The resulting xLm an d fit parameters are given in the last three 
columns of Table [l] for the three source models. The value of xLm now appears to be 
too small primarily due to the fact that 97-BLG-28 is less crowded than a typical star 
of its brightness, so that the frame quality undercertainties are overestimates of its true 
photometric scatter. What is important to note is that the uniformly bright source is again 
strongly ruled out and that all model parameters are left almost entirely unaffected by this 
alternate selection and treatment of the data, demonstrating the model's robustness. 

Both of the simpler models are special cases of the LD2 family of models: the LD1 
model requires the limb-darkening parameter d\ = in both / and V, while the UB model 
sets c\ = d\ = in both bands. This means that our LD1 and UB models, which correspond 
to points on the \ 2 hypersurface over restricted portions of the LD2 parameter space, 
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also correspond to points on the full LD2 x 2 hypersurface. The probability of obtaining 
a renormalized x 2 that is larger than that of the LD1 solution is > 10% for both the 
DoPHOT and frame quality methods of estimating relative photometric uncertainties; the 
probability of obtaining a renormalized x 2 larger than that of the UB is negligibly small. 

We conclude that the 2-parameter limb- darkening model is clearly superior to the 
uniform source model and marginally superior to the 1-parameter limb-darkening model. 
We now adopt it as our best model and use it to derive the physical parameters for the 
lens and the surface brightness profile of the source, which we now compare to expectations 
from stellar atmosphere models. 



4.2. Comparison to Stellar Atmospheric Models 

As we discussed in §0, both spectroscopic and photometric considerations indicate 
that the source star of 97-BLG-28 is a K giant in the Galactic bulge, with K2III being 
the most probable typing. The fits to our photometric data discussed in §|3] yield the 
limb-darkening coefficients c\ and d\ for the V and I bands separately. The corresponding 
surface brightness profiles for the source star, normalized so as to give a total flux equal 
to unity, are shown in Fig. 8 for our LD2 model of the high-quality photometric dataset 
using DoPHOT estimates for the relative photometric uncertainties. The surface brightness 
profiles derived from modeling the full dataset with empirical "frame quality" error bars 
produces nearly identical profiles (see coefficients in Table 1). Our results clearly indicate 
that, as expected, this giant is more limb-darkened in V than in /. 

Also shown in Fig. 8 are predictions from stellar atmospheric models for K0 and K5 
giants from two different groups flvanHamme 19931 piaz-Cordoves, Claret fc Gimenez 1995 



|Claret, Diaz-Cordoves fc Gimenez 1995|) both using the square-root limb-darkening law of 
Eq. 2. We have interpolated between the values given by these authors, labeling as K0 a 
star with effective temperature T = 4750 K and surface gravity log g = 2.15, and as K5 
a giant T = 4000 K and surface gravity log g = 1.75, in good agreement with standard 
convention ( [Lang 1992] ). The cooler, later type K5 giants are more limb-darkened than their 
warmer K0 cousins. Although the theoretical / profiles of both groups are nearly identical, 
the V profiles of Diaz-Cordoves, Claret & Gimenez (1995) are slightly more limb-darkened 
than those of van Hamme (1993). 

The surface brightness profiles derived from our photometric data modeled with 
the limb-darkening law of Eq. 2 are in excellent agreement with these predictions from 
atmospheric models. Although Fig. 8 indicates that the source star of 97-BLG-28 may be 
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slightly more limb-darkened in the I band than current models predict for its spectral type, 
this conclusion should be treated with considerable caution given the uncertainties in the 
light curve modeling, spectral typing, and the rather ad hoc nature of the fitting function 
for the surface brightness profile. 

Our results for 97-BLG-28 represent the first time that limb-darkening coefficients have 
been measured for a microlensing source. Limb-darkening was indicated but not measured 
by MACHO observations of the high-amplitude single-lens event MACHO 95-BLG-30; 
fixing limb-darkening coefficients to agree with stellar atmosphere models did produce a fit 
to the data that was marginally better (normalized A% 2 = 3.4) than uniform disk models 
(|Alcock et al. 1997b|). Our limb-darkening measurement for 97-BLG-28 is among only 



a very few determinations for normal giants, and respresents the first surface brightness 
profile for any star as distant as the Galactic bulge. The ability of PLANET photometry 
to yield measurements of limb-darkening coefficients in MACHO 97-BLG-28 is due to 
the overall characterization of the complete light curve combined with excellent temporal 
coverage (sampling times <, 30 minutes) over peak and limb caustic egress regions. 



5. Physical Parameters of the Lens 



In caustic crossing events, dense photometric monitoring allows not only the 
determination of t^, the time required for the lens to travel an angular distance 9e, but 
also a second time scale: the time £* required to travel the angular source radius For 
uniform rectilinear motion, the ratio £*/iE directly yields the angular size of the source 
in units of the Einstein ring, = 9*/9#. If 9* can be estimated by other means (e.g., 
photometric or spectroscopic typing), both 9e and the relative lens-source proper motion 
fi = #e/£e can thus be determined, yielding important information on lens kinematics 
flGould 1994| [Nemiroft fc Wickramasinghe 1994| |Witt fc Mao 19331 ; |Peng 1997| ). To date, 
proper motions derived from this technique have been published for only two binary events 
in the Magellanic Clouds, MACHO LMC-09 flAlcock et al. 1997a| ; [Bennett et al. 1995|) and 
MACHO 98-SMC-01 (|Afonso et al. 1998j |Albrow et al. 1999| ; Alcock et al.^999) ; |Udalski et 
al. 1998| ; |Rhie et al. T999|) , and one single-lens (point-caustic) event MACHO 95-BLG-30 
in the Galactic bulge ( [Alcock et al. 1997E| ). 



The relative proper motion can be written in terms of the model parameters and p*, 
and the physical source size R* and distance D$ as: 

R 



which for D$ = 8.0 kpc and R* 



P* tE 

15 Rq, as indicated by the spectral typing of §2.3, yield 
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a relative proper motion 

"= 19 ' 4 fe) "--^-' = 11.2 (^) .asday-. (4) 
This corresponds to a perpendicular velocity in the lens plane 

v = x D s /i = 155 a; ( — ) kms -1 , (5) 

V 15 R oJ 

where x = D L /D S . The (total) mass of the lens is given by 

so that for a disk lens halfway to the Galactic Center (x = 0.5), one obtains v — 78 kms -1 
and total lens mass M = 0.09 Mq, whereas a lens embedded in the Bulge with x = 0.9 
yields v = 140 kms" 1 and M = 0.81 Mq. Either value is consistent with Galactic 
kinematics. 

Assuming a source distance Ds = 8kpc and source radius in the range 13 — 17 R & 
(§2.3), the top panel of Fig. 9 gives the mass of the individual lensing components as a 
function of x. Since the binary mass ratio is q = 0.234, it is quite likely that the lighter of 
the lens objects has the mass of an M dwarf or less; if the lens lies in the Galactic disk with 
x < 0.8, at least the smaller of the binary components is a brown dwarf. If, on the other 
hand, the lensing binary resides in the bulge x > 0.8, the masses of its components would 
be consistent with that of a typical lower main-sequence binary with a projected separation 
a p = R*xd/ 'p* between 1 and 2 AU (bottom panel, Fig. 9). Disk lenses would have smaller 
separations. 



5.1. Blended Light: Light from the Lens Itself? 

Our model fits to 97-BLG-28 yield a fraction 1 — / of blended (non-lensed) light that 
is larger in the V than in the / band, indicating that the blend star is very much bluer 
than the background source star (Table 1). If the blend is due to a single star, it lies in 
an underpopulated portion of the CMD (Fig. 3), either because it is not at the mean field 
distance and has been improperly dereddened, or because it is in a short-lived phase of its 
evolution. A main sequence star with (V — 1) ^ 0.6 has I < 3.7, and thus is so intrinsically 
bright that it cannot appear as the blend anywhere along the reddening vector of Fig. 3 
unless it lies substantially behind the bulge, where it is likely to have more reddening than 
the 97-BLG-28 field, not less. 



-16- 



If, on the other hand, the blend star is somewhat in front of the bulge and experiences 
slightly less extinction, it could lie on the portion of the isochrone at (V — I) = —0.4 that 
corresponds to the transition zone between planetary nebulae and white dwarfs. Such an 
interpretation would be consistent, furthermore, with the blended light coming from the 
lens itself: a white dwarf mass of Mwd ~ 0.55 M© is consistent with the mass that the 
larger lens component (Mi) would have at a distance x ~ 0.88, slightly in front of the 
bulge. This dense portion of the Milky Way is a likely region for microlenses and should 
have an extinction similar to (but slightly less than) that of the field mean (see Fig. 3). 
If this explanation is correct, then the smaller lensing component would have a mass just 
above the hydrogen-burning limit at M 2 ~ 0.13 M . These conclusions are valid only if the 
blended light sensed by the light curve fitting is due to a single star, but any explanation 
requiring two stars to produce the color of the blend would be even more ad hoc. 

6. Conclusions 

The temporal coverage of the bulge microlensing event MACHO 97-BLG-28 from three 
PLANET observing sites longitudinally distributed in the southern hemisphere has resulted 
in an excellent characterization of the light curve in the / and V bands. The full data 
set consists of 513 photometric measurements in the /-band and 183 measurements in V, 
with 586 of these being of very high quality. During a 30-hour period over the peak of 
the light curve, sampling times between ~ 3 and 30 minutes were continuously maintained 
from the three sites. Less dense photometric monitoring over other portions of the light 
curve, combined with baseline measurements taken nearly a year later have resulted in 
the best-characterized microlensing light curve to date. Our spectrum of 97-BLG-28, 
taken at moderate magnification, indicates that the source star is likely to be a highly 
reddened K2 giant of solar metallicity. This conclusion is supported by a comparison of the 
color-magnitude diagram of the field to theoretical isochrones. 

Modeling of the PLANET photometric dataset for 97-BLG-28 clearly indicates that 
the disk of the resolved source is transited by the cusp of the central caustic generated by a 
binary lens with mass ratio q = 0.23 and instantaneous separation, in units of the Einstein 
ring radius, of d = 0.69. A source of uniform surface brightness is strongly ruled out by 
our data; the derived square root limb-darkening coefficients in the I and V bands derived 
from the photometric data alone are in excellent agreement with expectations from stellar 
atmospheric models for K giants. Under the assumption that the source star is a K2III, the 
models appear slightly more limb-darkened in I than our measurements, but the difference 
may not be significant. 
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All of the conclusions above are independent of the (unknown) lens-source distance 
ratio x and (rather well-constrained) source distance D s and radius R*. Assuming a source 
in the Galactic bulge at 8kpc of radius R* = 15 ± 2 R Q , as indicated by our spectral typing, 
the total mass of the binary and its projected separation can be determined. If the lens 
resides in the bulge, it is likely to be a lower main sequence binary with both components 
above the hydrogen burning limit and a projected separation between 1 and 2 AU. If the 
lens resides in the disk with the lens-source distance ratio x < 0.8, the projected separation 
would be smaller and one or both of the components would have a mass in the brown dwarf 
regime. Assuming the same source distance and radius, the relative lens-source proper 
motion derived from our modeling is /i = 19.4 ± 2.6 kms -1 kpc -1 or fi = 11.2 ± 1.5 /ias/day, 
consistent with the lens having disk or bulge kinematics. The uncertainty in the proper 
motion is dominated by the uncertainty in the source radius, not by the characterization 
of the microlensing light curve. If the unlensed, blended light indicated by our models is 
due to a single star, its color and magnitude suggest that it is a young white dwarf in the 
bulge. This conclusion is consistent with the light coming from the lens itself: the larger 
lens component would have a typical white dwarf mass at a bulge distance of x ~ 0.88, in 
which case its companion would be a very faint late M dwarf. 

These results mark the first time that the surface structure of a source star has been 
measured via microlensing and the first limb-darkening determination of a star as distant 
as the Galactic Bulge by any technique. As the number of caustic-crossing events alerted 
in real time continues to rise, and the quality and frequency of photometric monitoring of 
microlensing light curves continues to improve, microlensing holds the promise of making 
a substantial contribution to the field of stellar atmospheres by determining the surface 
profiles of normal stars too faint or small to be easily measured by other techniques. 

PLANET thanks the MACHO collaboration for providing the original real-time 
electronic alert of this event; such alerts are crucial to the success of our intensive 
microlensing monitoring. We are grateful to observers from the Astronomical Society of 
Tasmania, especially Bob Coghlan, who made many of the key observations reported here. 
We thank Andy Gould for reading a penultimate version of this manuscript. Financial 
support from Nederlands Wetenschapelijk Onderzoek, through award ASTRON 781.76.018, 
is gratefully acknowledged. PLANET members also wish to thank The Leids Sterrewacht 
Foundation, the South African Astronomical Observatory, Canopus Observatory, the 
European Southern Observatory, and Perth Observatory for the generous allocations of 
time that make our work possible. The work of M.D. has been financed by a research grant 
from the Deutsche Forschungsgemeinschaft while at STScI, Baltimore and by a Marie Curie 
Fellowship at Kapteyn Institute, Groningen. 
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Table 1. PLANET model parameters for Galactic bulge event MACHO 97-BLG-28. 





Formal DoPHOT Error Bars 


"Frame Quality" 


Error Bars 




LD2 


LD1 


UB 


LD2 


LD1 


UB 


V 2 . 

A mm 


1 01 Q 

lyio 


lyoU 


ozOO 


OP A 
Z04 


Ofi 1 
ZOl 


450 


# d.o.f 


00 ( 


kpo 

ooy 


C71 
OIL 


PP1 
00 ( 


PPO 

ooy 


671 


# parameters 


1 o 


1 7 


10 


1 o 

iy 


1 7 


15 


*e [d] 


07 Q 
Z / .0 


07 o 

z / .z 


OP 1 
ZO. 1 


Z 1 .0 


07 Q 

At .6 


26.1 


to [d] 


oner pro 

oyo.oo 


oyo.oo 


one po 

oyo.oo 


oyo.oo 


one; kq 

oyo.oo 


895.63 




U.UUzy 


U.UUoU 


n noi 
U.Uzl 


n nnon 

u.uuzy 


n nriQ/i 
U.UUo4 


0.021 


d 


U.ooo 


U.00 ( 


n £?7o 
U.o to 


n esc 
U.ooo 


n coo 
U.ooo 


0.680 




U.Z04 


n oqo 
U.ZoZ 


n 077 
U.Z / i 


U.Z04 


n oqi 
U.zOl 


0.276 


ct 


1 /I OP 
1.4Z0 


1 AO A 

1.4Z4 


1.4U0 


1 A 07 
1.4Z / 


1 A OO 
1.4ZZ 


1.419 


/i 


U.y I 


n nc 

u.yo 


n oq 

u.yo 


n 07 

u.y ( 


n no 

u.yo 


0.96 


fy 


0.84 


0.84 


0.85 


0.85 


0.88 


0.87 


P* 


0.0286 


0.0287 


0.0296 


0.0286 


0.0289 


0.0295 




0.40 


0.75 





0.38 


0.83 





cv 


0.55 


0.92 





0.55 


0.95 





dj 


0.37 








0.38 








dy 


0.44 








0.43 










15.632 


15.633 


15.620 


15.630 


15.631 


15.597 




15.643 


15.644 


15.631 


15.643 


15.645 


15.609 


TI 


15.708 


15.703 


15.727 


15.697 


15.697 


15.673 


mg V 


17.899 


17.898 


17.889 


17.909 


17.930 


17.900 


sv 


17.906 


17.905 


17.903 


17.914 


17.938 


17.903 


TV 

m lv 


17.926 


17.918 


17.965 


17.947 


17.959 


17.968 



Note. In separate columns, model parameters for 2-parameter limb- 

darkened (LD2), 1-parameter limb-darkened (LD1) and uniformly bright (UB) 
source models are listed. The first three columns indicate fits of high quality data 
with DoPHOT-reported error bars; the last three fits are of all data with empirical 
"Frame Quality" uncertainties based on typical scatter for similar brightness stars 
on individual frames. The baselines mo refer to the total modeled flux (source + 
blend) at the given observing site (L=LaSilla, S=SAAO, T=Tasmania) and in the 
given band (I or V). 
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Figure Captions 

Fig. 1. — Left: The high-quality subset of PLANET photometry for 97-BLG-28 spanning a 
300 day period and consisting of 431 /-band points (bottom) and 155 V-band points (top) 
from three PLANET observing stations in Tasmania (green), Chile (black), and South Africa 
(red) is shown with baseline points taken the following season. Light curves for our best- 
fitting limb-darkened model (LD2) are superimposed. Right: An enlargement showing the 
30-day period near the peak. 

Fig. 2. — Flux-calibrated, combined red and blue spectra taken with EFOSC on the ESO 
3.6m. The spectrum has not been dereddened, and was taken before the anomalous peak 
when the event was only slightly magnified. 

Fig. 3. — Calibrated, dereddened color-magnitude diagram of the MACHO 97-BLG-28 field, 
taken from SAAO lm PLANET data. A reddening of E(V — I) = 1.095 has been applied, as 
spectral typing and isochrone fitting to the red clump and red giant branch. Also shown are 
the dereddened position on the diagram of the unmagnified, blended source (filled circle), 
and the positions of the blend and unblended source (open circles) as inferred from modeling. 
Isochrones from Bertelli et al. 1994 of age 6 Gyr and solar and super-solar metallicity are 
overplotted. The faint blend star would lie along indicated reddening vector (i.e., would be 
fainter and redder than shown) if it suffered less extinction than the field mean. 

Fig. 4. — Left: Light curves for our best fit, the 2-parameter limb-darkening model (LD2), 
superimposed on PLANET V and / data sets for a 3-day period centered on the cusp- 
crossing (bottom) and a 16-hour period during which the stellar limb swept over the cusp 
(top). Right: Same for the uniform bright source model (UB). DoPHOT estimates for 
relative photometric uncertainties have been assumed; the rescaled error bars are generally 
comparable to the size of the plotted points, and so have been omitted. 

Fig. 5. — Configuration of binary lens and trajectory of the extended source on the sky from 
the 2-parameter limb-darkening model. The lens caustics are shown as the closed, pointed 
shapes in the center and right side of the diagram. The path of the source is indicated by 
the straight parallel lines; the distance between the two outer parallel lines corresponds to 
the source width. Lens positions are shown as crosses; Mi is the heavier of the two. All 
distances are in units of the angular Einstein radius. 

Fig. 6. — Top: Magnification curves for 10 uniformly bright, equal-area rings of a source 
crossing a cusp with the geometry of Fig. 5. Bottom: The contribution to the total light at 
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any given time from each of the uniformly bright rings. 

Fig. 7. — V and I band residuals (Model - Data) for the models of Fig. 4 during the brief 
period in which the cusp swept over the trailing limb of the source. (Rescaled) formal 
DoPHOT error bars are shown. The different sampling rates and observing conditions at 
the three sites (successively La Silla, Tasmania and SAAO) is apparent. Left: Residuals for 
our best fitting 2-parameter limb-darkening model (LD2). Right: Same for the uniformly 
bright source model (UB). 

Fig. 8. — Source radial profiles in the I-band (left) and V-band (right) for uniformly-bright 
(thin solid line) and our 2-parameter limb-darkened source (thick solid line) models. Also 
shown are the theoretical expectations in the Cousins I and Johnson V bands based on 
atmospheric models for KO — K5 giants (dashed lines) from two different literature sources 
(see text). 

Fig. 9. — Top panel: The larger (Mi, solid lines) and smaller (M 2 , dashed lines) lens masses 
shown as a function of x, the ratio of lens to source distances. Parallel lines enclose the range 
of uncertainty dominated by the uncertainty in the source radius, estimated to lie between 
13 and 17 Rq. A source distance of 8kpc has been assumed. Bottom panel: The projected 
binary separation is shown in AU as a function of x for the same source radii. 
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